Sideband Cooling Micromechanical Motion to the Quantum Ground State 
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The advent of laser cooling techniques revolutionized 
the study of many atomic-scale systems. This has fu- 
eled progress towards quantum computers by prepar- 
ing trapped ions in their motional ground state Q]> 
and generating new states of matter by achieving Bose- 
Einstein condensation of atomic vapors [2|. Analogous 
cooling techniques 1 3 , 4 1 provide a general and flexible 
method for preparing macroscopic objects in their mo- 
tional ground state, bringing the powerful technology of 
micromechanics into the quantum regime. Cavity opto- 
or electro-mechanical systems achieve sideband cooling 
through the strong interaction between light and mo- 
tion Il5l-[l5l. However, entering the quantum regime, 
less than a single quantum of motion, has been elu- 
sive because sideband cooling has not sufficiently over- 
whelmed the coupling of mechanical systems to their 
hot environments. Here, we demonstrate sideband cool- 
ing of the motion of a micromechanical oscillator to the 
quantum ground state. Entering the quantum regime 
requires a large electromechanical interaction, which is 
achieved by embedding a micromechanical membrane 
into a superconducting microwave resonant circuit. In 
order to verify the cooling of the membrane motion 
into the quantum regime, we perform a near quantum- 
limited measurement of the microwave field, resolving 
this motion a factor of 5.1 from the Heisenberg limit 
|3|. Furthermore, our device exhibits strong-coupling 
allowing coherent exchange of microwave photons and 
mechanical phonons lHol . Simultaneously achieving 
strong coupling, ground state preparation and efficient 
measurement sets the stage for rapid advances in the 
control and detection of non-classical states of motion 
O 2HL possibly even testing quantum theory itself in 
the unexplored region of larger size and mass [19|. The 
universal ability to connect disparate physical systems 
through mechanical motion naturally leads towards fu- 
ture methods for engineering the coherent transfer of 
quantum information with widely different forms of 
quanta. 

Mechanical oscillators that are both decoupled from 
their environment (high quality factor Q) and placed in the 
quantum regime could allow us to explore quantum me- 
chanics in entirely new ways lfT7H2Tl . For an oscillator to 
be in the quantum regime, it must be possible to prepare 
it in its ground state, to arbitrarily manipulate its quantum 
state, and to detect its state near the Heisenberg limit. In or- 
der to prepare an oscillator in its ground state, its tempera- 



ture T must be reduced such that k B T < h£l m , where Q m 
is the resonance frequency of the oscillator, k B is Boltz- 
mann's constant, and h is the reduced Planck's constant. 
While higher resonance frequency modes (> 1 GHz) can 
meet this cooling requirement with conventional refriger- 
ation (T < 50 mK), these stiff oscillators are difficult to 
control and to detect within their short mechanical life- 
times. One unique approach using passive cooling has suc- 
cessfully overcome these difficulties by using a piezoelec- 
tric dilatation oscillator coupled to a superconducting qubit 
J22). Unfortunately, this method is incompatible with the 
broad range of lower frequency, high Q, flexural mechani- 
cal modes. In order to take advantage of the attractive me- 
chanical properties of these oscillators, an alternative active 
cooling method is required, one that can reduce the oscil- 
lator's temperature below that of the surrounding environ- 
ment. 

Cavity opto- or electro-mechanical systems [4] natu- 
rally offer a method for both detecting mechanical mo- 
tion and cooling a mechanical mode to its ground state 
ll23l l24l . An object whose motion alters the resonance 
frequency u c of an electromagnetic cavity experiences a 
radiation pressure force governed by the parametric inter- 
action Hamiltonian: H mt = hGhx, where G = du c /dx, 
h is the cavity photon number, and x is the displacement 
of the mechanical oscillator. By driving the cavity at a 
frequency u)<±, the oscillator's motion produces upper and 
lower sidebands at Wd db f2 m . Because these sideband pho- 
tons are inelastically scattered from the drive field, they 
provide a way to exchange energy with the oscillator. If 
the drive field is optimally detuned below the cavity res- 
onance A = u d — lo c = — f2 m , photons will be prefer- 
entially up-converted to lo c because the photon density of 
states is maximal there (Fig lb). When an up-converted 
photon leaves the cavity, it removes the energy of one me- 
chanical quantum (one phonon) from the motion. Thus, 
the mechanical oscillator is damped and cooled via this 
radiation-pressure force. Because the mechanical motion is 
encoded in scattered photons exiting the cavity, a quantum- 
limited measurement of this photon field provides a near 
Heisenberg-limited detection of mechanical motion J25j. 

While there has been substantial progress in cooling me- 
chanical oscillators with radiation pressure forces, side- 
band cooling to the quantum mechanical ground state has 
been an outstanding challenge. Cavity optomechanical sys- 
tems have realized very large sideband cooling rates J8]- 
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Q2KT41Q3); however, these rates are not sufficient to over- 
come the larger thermal heating rates of the mechanical 
modes. Because electromechanical experiments use much 
lower-energy photons 151(71 [l"3ll , they are naturally compat- 
ible with operation below 100 mK, but have consequently 
suffered from weak electromechanical interactions and in- 
efficient detection of the photon fields. 

Here, we present a cavity electromechanical system 
where a flexural mode of a thin aluminum membrane is 
parametrically coupled to a superconducting microwave 
resonant circuit. Unlike previous microwave systems, this 
device achieves large electromechanical coupling by con- 
centrating nearly all the microwave electric fields near the 
mechanical oscillator |[l"6l . The oscillator is a 100 nm thick 
aluminum membrane with a diameter of 15 um, suspended 
50 nm above a second aluminum layer on a sapphire sub- 
strate [261 ( see Fig. 1). These two metal layers form a 
variable parallel-plate capacitor that is shunted by a 12 nH 
spiral inductor. This combination of capacitor and induc- 
tor creates a microwave cavity whose resonance frequency 
depends on the mechanical displacement of the membrane 
and is centered at uj c = 2ir x 7.54 GHz. The device is oper- 
ated in a dilution refrigerator at 15 mK, where aluminum is 
superconducting, and the microwave cavity has a total en- 
ergy decay rate of k ~ 2ir x 200 kHz. As expected from 
the dimensions of the membrane, tl m = 2ir x 10.56 MHz, 
and we find an intrinsic damping rate of T m = 2n x 32 Hz, 
resulting in a mechanical quality factor Q m = Jl m /T m = 
3.3 x 10 5 . The oscillator mass m = 48 pg implies that the 
zero-point motion is x zp = y/h/(2mn m ) = 4.1 fm. With 
a ratio of £l m / K > 50, our system is deep in the resolved- 
sideband regime and perfectly suited for sideband cooling 
to the mechanical ground state Il23ll24l . 

To measure the mechanical displacement, we apply a mi- 
crowave field, which is detuned below the cavity resonance 
frequency by A = — fi m > through heavily attenuated coax- 
ial lines to the feed line of our device. The upper sideband 
at oj c is amplified with a custom-built Josephson parametric 
amplifier (JPA) [27, 28 1 followed by a low-noise cryogenic 
amplifier, demodulated at room temperature, and finally 
monitored with a spectrum analyzer. The thermal motion 
of the membrane creates an easily resolvable peak in the 
microwave noise spectrum. As described previously [ 28 1, 
this measurement scheme constitutes a nearly shot-noise- 
limited microwave interferometer with which we can mea- 
sure mechanical displacement with minimum added noise 
close to fundamental limits. 

In order to calibrate the demodulated signal to the mem- 
brane's motion, we measure the thermal noise spectrum 
while varying the cryostat temperature (Fig. lc). Here a 
weak microwave drive (~ 3 photons in the cavity) is used 
in order to ensure that radiation pressure damping and cool- 
ing effects are negligible. When O m 3> K S> T m and 
A = — Q m , the displacement spectral density S x is re- 
lated to the observed microwave noise spectral density S 
by: S x = 2(K£l m /GK cx ) 2 S '/ P Q , where k cx is the cou- 



pling rate between the cavity and the feed line, and P a is 
the power of the microwave drive at the output of the cav- 
ity. According to equipartition, the area under the reso- 
nance curve of displacement spectral density S x must be 
proportional to the effective temperature of the mechani- 
cal mode. This calibration procedure allows us to convert 
the sideband in the microwave power spectral density to 
a displacement spectral density and to extract the thermal 
occupation of the mechanical mode. In Fig. lc we show 
the number of thermal quanta in the mechanical resonator 
as a function of T. The linear dependence of the inte- 
grated power spectral density with temperature shows that 
the mechanical mode equilibrates with the cryostat even 
for the lowest achievable temperature of 15 mK. This tem- 
perature corresponds to a thermal occupancy n m = 30, 
where n m = [exp(H{l m /k B T) — The calibra- 

tion determines the electromechanical coupling strength 
G/2tt = 49 ± 2 MHz/nm. With the device parameters, 
we can investigate both the fundamental sensitivity of our 
measurement as well as the effects of radiation pressure 
cooling. 

The total measured displacement noise results from two 
sources: the membrane's actual mean-square motion S x h 
and the apparent motion S x mp due to imprecision of the 
measurement. Fig. 2a demonstrates how the use of low- 
noise parametric amplification significantly lowers S x mp , 
resulting in a reduction in the white-noise background by a 
factor of more than 30. This greatly increases the signal-to- 
noise ratio of the membrane's thermal motion, reducing the 
required integration time to resolve the thermal peak by a 
factor of 1000. To investigate the measurement sensitivity 
in the presence of dynamical backaction, we regulate the 
cryostat temperature at 20 mK and increase the amplitude 
of the detuned microwave drive while observing modifica- 
tions in the displacement spectral density. We quantify the 
strength of the drive by the resulting number of photons n d 
in the microwave cavity. As shown in Fig. 2b, the measure- 
ment imprecision 5*™ p is inversely proportional to n&. At 
the highest drive power (rid ~ 10 5 ), the absolute displace- 
ment sensitivity is 5.5 x 10~ 34 m 2 /Hz. 

As expected, the increased drive power also damps and 
cools the mechanical oscillator Il3l l23ll24l . The total me- 
chanical dissipation rate T' m is the sum of the intrinsic dis- 
sipation r m and the radiation-pressure -induced damping 
resulting from scattering photons to the upper/lower side- 
band r = r + - r_, where T± = 4 9 2 k/[k 2 + 4(A ± 
O m ) 2 ]. Here, g is the coupling rate between the cavity 
and the mechanical mode, which depends on the ampli- 
tude of the drive: g = Gx zp ^/rTd. Fig. 2c shows the 
measured values of k, g and T' m as the drive increases. 
The radiation-pressure damping of the mechanical oscil- 
lator becomes pronounced above a cavity drive amplitude 
of approximately 75 photons, at which point T = T m and 
the mechanical linewidth has doubled. 

While the absolute value of the displacement impre- 
cision decreases with increasing power, the visibility of 
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the thermal mechanical peak no longer improves once 
the radiation-pressure force becomes the dominant dissi- 
pation mechanism for the membrane. By expressing the 
imprecision as equivalent thermal quanta of the oscillator 
"-imp = F' m S'% ap /8xl p , we see that the visibility of the 
thermal noise above the imprecision no longer improves 
once the drive is much greater than n d ~ 100 (Fig. 2d). 
This is because a linear decrease in S l ™ p is balanced by 
a linear increase in T' m due to radiation-pressure damp- 
ing. The asymptotic value of ra imp is a direct measure of 
the efficiency of the microwave measurement. Ideally, for 
a lossless circuit, a quantum-limited microwave measure- 
ment would imply n imp = 1/4. The incorporation of the 
low-noise JPA improves n imp close to this ideal limit, re- 
ducing the asymptotic value of re imp from 70 to 1.9 quanta. 
This level of sensitivity is crucial, as we will now use this 
measurement to resolve the residual thermal motion of the 
membrane as it is cooled into the quantum regime. 

Beginning from a cryostat temperature of 20 mK and a 
thermal occupation of = 40 quanta, the fundamen- 
tal mechanical mode of the membrane is cooled by the 
radiation-pressure forces. Figure 3a shows the displace- 
ment spectral density of the motional sideband as n d is 
increased from 18 to 4,500 photons along with fits to a 
Lorentzian lineshape (shaded area). As described above, 
this increased drive results in three effects on the spectra: 
lower noise floor, wider resonances and smaller area. As it 
is the area that corresponds to the mean-square motion of 
the membrane, it directly measures the effective tempera- 
ture of the mode. At a drive intensity that corresponds to 
4,000 photons in the cavity, the thermal occupation is re- 
duced below one quantum of mechanical motion, entering 
the quantum regime. 

Observing the noise spectrum over a broader frequency 
range reveals that there is also a second Lorentzian peak 
with linewidth k whose area corresponds to the finite ther- 
mal occupation n c of the cavity. Over a broad frequency 
range it is no longer valid to evaluate the cavity parame- 
ters at a single frequency to infer the spectrum in units of 
S x . Instead, Fig. 3b shows the noise spectrum in units of 
sideband power normalized by the power at the drive fre- 
quency, S/ P Q . These two sources of noise originating from 
either the mechanical or the electrical mode interfere with 
each other and result in noise squashing |[T3l and eventu- 
ally normal-mode splitting 11291 once 2g > Us- 
ing a quantum-mechanical description applied to our cir- 
cuit |fT3ll25l , the expected noise spectrum is 

1 2k ox [«n c (r^ + 4(5 2 ) + 4r m n^g 2 ] 
S/frw = - + n add + — -g 

|V+(« + 2i(* + A)J(r m + 2i*)| 

(!) 

where 5 = uj — O m , A = u d + Q m — u c , and ra add 
is added noise of the microwave measurement expressed 
as an equivalent number of microwave photons. Fig. 3b 
shows the measured spectra and corresponding fits (shaded 
region) to Eq. 1 as the electromechanical system evolves 



first into the quantum regime (n m , n c < 1) and then into 
the strong-coupling regime (2g > ft/2). The results are 
summarized in Fig. 3c, where the thermal occupancy of 
both the mechanical and electrical modes are shown as a 
function of re d . For low drive power, the cavity shows no 
resolvable thermal population (to within our measurement 
uncertainty of 0.05 quanta) as expected for a 7.5 GHz mode 
at 20 mK. While it is unclear whether the observed popula- 
tion at higher drive power is a consequence of direct heat- 
ing of the substrate, heating of the microwave attenuators 
preceding the circuit, or intrinsic cavity frequency noise, 
we have determined that it is not the result of frequency or 
amplitude noise of our microwave generator, as this noise 
is reduced far below the microwave shot-noise level with 
narrow-band filtering and cryogenic attenuation (see Sup- 
plementary Information). Sideband cooling can never re- 
duce the occupancy of the mechanical mode below that of 
the cavity. Therefore, in order for the system to access the 
quantum regime, the thermal population of the cavity must 
remain less than one quantum. Assuming f2 m 3> n, the 
final occupancy of a mechanical mode is fl29l 

n =n T f^ 4g2 + K2 V n( ^ \ (2) 

This equation shows that for moderate coupling (\/nT m <C 
g <C k) the cooling of the mechanical mode is linear in the 
number of drive photons. Beyond this regime, the onset of 
normal-mode splitting abates further cooling. Here the me- 
chanical cooling rate becomes limited not by the coupling 
between the mechanical mode and the cavity, but instead 
by the coupling rate k between the cavity and its environ- 
ment [29]. Thus, the final occupancy of the mechanical 
mode can never be reduced to lower than n^r m / n, and 
a stronger parametric drive will only increase the rate at 
which the thermal excitations Rabi oscillate between the 
cavity and mechanical modes. For our device we achieve 
the desired hierarchy: as the coupling is increased, we first 
cool to the ground state and then enter the strong-coupling 
regime (n^T m < k < g). Once re d exceeds 2 x 10 4 , the 
mechanical occupancy converges toward the cavity popu- 
lation, reaching a minimum of 0.34 ± 0.05 quanta. At the 
highest power drive power (n d = 2 x 10 5 ) the mechan- 
ical mode has hybridized with the cavity, resulting in the 
normal-mode splitting characteristic of the strong-coupling 
regime |[T6l . This level of coupling is required to utilize 
the hybrid system for quantum information processing, as 
it is only in the strong-coupling regime that a quantum state 
may be manipulated faster than it decoheres from the cou- 
pling of either the electromagnetic or mechanical modes to 
the environment. 

Together the measurements shown in Fig. 2 and 3 quan- 
tify the overall measurement efficiency of the system. The 
Heisenberg limit requires that a continuous displacement 
measurement is necessarily accompanied by a backaction 

force J3j[12l|25j, such that \/ ST^Sp > ft, where Sp is the 
force noise spectral density. From the thermal occupancy 
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and damping rate of the mechanical mode, we extract the 
total force spectral density Sp = 4M2 m mr^ n (n in + 1/2). 
This places a conservative upper bound on the quantum 
backaction by assuming that it alone is responsible for 
the finite occupancy of the mechanical mode. This ex- 
periment achieves the closest approach to Heisenberg- 
limited displacement detection to date lfT4l l25ll with 
a lowest imprecision-backaction product 
4hy/n imp (n m + 1/2) = (5.1 ± 0A)h. Thus, this me- 
chanical device simultaneously demonstrates ground-state 
preparation, strong-coupling and near quantum-limited de- 
tection. 

Looking forward, this technology offers a feasible route 
to achieve many of the longstanding goals for quantum me- 
chanical systems. These prospects include a direct mea- 
surement of the zero-point motion, observation of the fun- 
damental asymmetry between the rate of emission and ab- 
sorption of phonons [ 1 1, quantum nondemolition measure- 
ments J3j and generation of entangled states of mechanical 
motion lfT71[T8l . Furthermore, combining this device with a 
single-photon source and detector (such as a superconduct- 
ing qubit li22l [30)) would enable preparation of arbitrary 
quantum states of mechanical motion as well as observa- 
tion of a single excitation as it Rabi oscillates between a 
7 GHz photon and a 10 MHz phonon EDI . Because the 
interaction between the mechanical mode and the cavity 
is parametric, the coupling strength is inherently tunable 
and can be turned on and off quickly. Thus, once a quan- 
tum state is transfered into the mechanical mode, it can be 
stored there for a time r t h = l/(n^L m ) > 100 us be- 
fore absorbing one thermal phonon from its environment. 
As this timescale is much longer than typical coherence 
times of superconducting qubits, mechanical modes offer 
the potential for delay and storage of quantum information. 
Lastly, because mechanical oscillators can couple to light 
of any frequency, they could serve as a unique intermediary 
that transfers quantum information between the microwave 
and optical domains ETTl . 

These measurements demonstrate the power of sideband 
techniques to cool a macroscopic (~ 10 12 atoms) mechan- 
ical mode, beyond what is feasible with conventional re- 
frigeration techniques, into the quantum regime. These 
broadly applicable methods for state preparation, manipu- 
lation and detection, pave the way to access the quantum 
nature of a wide class of long-lived mechanical oscilla- 
tors. Through the strong interaction between photons and 
phonons, mechanical systems can now inherit the experi- 
mental and theoretical power of quantum optics, opening 
the field of quantum acoustics. 
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Cryostat Temperature (mK) 

FIG. 1 . Schematic description of the experiment, a, Colorized scanning electron micrograph showing the aluminum (grey) elec- 
tromechanical circuit fabricated on a sapphire (blue) substrate, in which a 15 urn diameter membrane is lithographically suspended 
50 nm above a lower electrode. The membrane's motion modulates the capacitance, and hence, the resonance frequency of the super- 
conducting microwave circuit, b, A coherent microwave drive inductively coupled to the circuit acquires modulation sidebands due to 
the thermal motion of the membrane. The upper sideband is amplified with a nearly quantum-limited Josephson parametric amplifier 
within the cryostat. c, The microwave power in the upper sideband provides a direct measurement of the thermal occupancy of the me- 
chanical mode, which may be calibrated in situ by varying the temperature of the cryostat. The mechanical mode shows thermalization 
with the cryostat at all temperatures, yielding a minimum thermal occupancy of 30 mechanical quanta without employing sideband- 
cooling techniques. The inset illustrates the concept of sideband cooling. When the circuit is excited with a detuned microwave drive 
such that A = — f2 m , the narrow line shape of the electrical resonance ensures that photons are preferentially scattered to higher energy, 
providing a cooling mechanism for the membrane. 
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FIG. 2. Displacement sensitivity in the presence of radiation-pressure damping, a, The displacement spectral density measured 
with (red) and without (blue) the Josephson parametric amplifier. As the parametric amplifier greatly reduces the total noise of the 
microwave measurement, the time required to resolve the thermal motion is reduced by a factor of 1000. b, As the microwave drive 
power is increased, the absolute displacement sensitivity, S 1 ™ 13 improves, reaching a minimum of 5.5 x 10~ 34 m 2 /Hz at the highest 
power, c, The parametric coupling g between the microwave cavity and the mechanical mode increases as y/n^. This coupling damps 
the mechanical mode from its intrinsic linewidth of r m = 2ty x 32 Hz until it is increased to that of the microwave cavity k. d, The 
relative measurement imprecision, in units of mechanical quanta, depends on the product of S™ 13 and T' m . Thus, once the power 
is large enough that radiation-pressure damping overwhelms the intrinsic mechanical dissipation, n lmp asymptotically approaches a 
constant value (n imp — 1.9), which is a direct measure of the overall efficiency of the photon measurement. 
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FIG. 3. Sideband cooling the mechanical mode to the ground state, a, The displacement noise spectra and Lorentzian fits (shaded 
region) for five different drive powers. With higher power, the mechanical mode is both damped (larger linewidth) and cooled (smaller 
area) by the radiation pressure forces, b, Over a broader frequency span, the normalized sideband noise spectra clearly show both the 
narrow mechanical peak and a broader cavity peak due to finite occupancy of the mechanical and electrical modes, respectively. A 
small, but resolvable, thermal population of the cavity appears as the drive power increases, setting the limit for the final occupancy of 
the coupled optomechanical system. At the highest drive power, the coupling rate between the mechanical oscillator and the microwave 
cavity exceeds the intrinsic dissipation of either mode, and the system hybridizes into optomechanical normal modes, c, Starting in 
thermal equilibrium with the cryostat at T = 20 mK, sideband cooling reduces the thermal occupancy of the mechanical mode from 
n m = 40 into the quantum regime, reaching a minimum of n m = 0.34 ± 0.05. These data demonstrate that the parametric interaction 
between photons and phonons can initialize the strongly coupled, electromechanical system in its quantum ground state. 
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Supplementary Information for "Sideband Cooling Micromechanical Motion to the Quantum 

Ground State" 



NOISE SPECTRUM OF AN OPTOMECHANICAL 
SYSTEM 

A mechanical degree of freedom that parametrically cou- 
ples to the cavity resonance frequency modifies the power 
emerging from the cavity by scattering photons to the up- 
per or lower mechanical sidebands. To calculate the full 
noise spectrum of the optomechanical system, we follow 
the general method of input-output theory (T). We define 
g = Gx zpy /n^, where G = duj c /dx, x zp = -^/h/2mfl m , 
m is the mass, oj c is the cavity resonance frequency, O m 
is the mechanical resonance frequency and rid is the num- 
ber of photons in the cavity due to a drive at frequency Wd- 
Furthermore, we define the response functions of the me- 
chanical and cavity modes as x^T 1 = k/2 + j(S + A) and 
X" 1 = r m /2 + jS, where r m is the mechanical dissipa- 
tion rate, k is the cavity dissipation rate, 5 = u — S7 m , 
A = w d — w c + and j = y/—l. K is total cavity dissi- 
pation rate due to both the intentional coupling to the trans- 
mission line k cx and the intrinsic losses Kq. From these pa- 
rameters, we define the optomechanical self-energy J2j [3j 
as a function of 5: 



m = ~39 2 [XM 

~ -jg 2 Xc{5) 



(SI) 
(S2) 



The approximation assumes that the drive is near the opti- 
mal detuning for cooling (|A| <C O m ) and the system is 
sufficiently in the good-cavity limit (S7 m 3> n) such that 
the cavity response at {5 + 2f2 m ) may be neglected. Now 
the effective mechanical response function x m including 
the optomechanical effects is: 



Xn 



1 + JXmS 

,-1 



a 2 + Xn^Xc 1 



(S3) 
(S4) 



The noise at the output of the cavity is characterized by 
the noise operator 6 out , which is related to the cavity field 
operator a by 6 out = f3K cx a. (3 is a dimensionless factor 
that depends on the geometry. Our circuit (shown schemat- 
ically in Fig. SI) couples power from the cavity equally 
to the output and back to the input so here j3 = 1/2. In 
principle, this fraction could be engineered by coupling 
asymmetrically to the input and the output, or by using a 
single port cavity (f3 = 1). Following directly the the- 
oretical analysis of previous work (5J 0), we consider the 
noise operators f) m and fj c associated with the mechanical 
and cavity modes respectively, which satisfy the relations 



( input Q ( output Q 

I K ex 



cavity 

CO. 



dissipation to environment 



FIG. S 1 . Cavity coupling block diagram. 



(VmVm) = n m an d (vtVc) = n c- Thus, the output noise is 

a 



6out = - VP^cxXcVk (1 - g 2 XmXc) fj c 

- VPKcxXcyfi\^(jgx m ) ^/m* 

In the frequency domain, the power spectral density of the 
noise at the output (in units of W/Hz) is 5 = hui(b out b on t), 

4^/3At cx (r^ + 45 2 )nn c 



S 



+ 



4g* + ( K + 2j(6 + A)) (T m + 2j5) 
16hoj(3K cx g 2 T m nl 1 



« + 2j(<5 + A) (r m + 2j,5) 



The first term simply represents the thermal noise of a 
cavity with occupancy n c whose spectral weight is dis- 
tributed over the 'dressed' cavity mode. The 'dressed' cav- 
ity mode includes the effect of optomechanically induced 
transparency J7J? ? ] and reduces to a single Lorentzian 
lineshape in the limit of weak coupling (g <C y/K,T m ). 
The second term is the thermal noise of the mechanical 
mode with its modified mechanical susceptibility. Unlike 
previous derivations [4|, we have not assumed the weak- 
coupling regime. Thus, as this equation is valid in both 
the weak- and strong-coupling regimes, it gives a unified 
description of the thermal noise spectrum even in the pres- 
ence of normal-mode splitting. Finally, the total noise at 
the output of the measurement including the vacuum noise 
of the photon field and the added noise of the measurement 
is 



S 1 , 

to = 2 +H » 



add 



+ 



4/3k cx [ K n c (r^ + 45 2 ) + 4T m nlg 2 } 



4gZ+( K + 2j(6 + A))(T m + 2j5) 

(S5) 
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FIG. S2. Detailed schematic diagram. A microwave generator creates a tone at the drive frequency. This signal is filtered with a 
resonant cavity at room temperature and split into two arms. The first arm excites the cavity through approximately 53 dB of cryogenic 
attenuation. In order to avoid saturating the low-noise amplifier with the microwave drive tone, the second arm is used to cancel 
the drive before amplification. A computer-controlled variable attenuator and phase shifter are run in a feedback loop to maintain 
cancellation at the part per million level. A second microwave generator is used to provide the pump tone for the Josephson parametric 
amplifier (JPA) as well as the reference oscillator for the mixer. This pump tone is 1.3 MHz above uj c so that the JPA is operated as 
a non-degenerate parametric amplifier, which measures both quadratures of the electromagnetic filed at the upper sideband frequency. 
The last stage of attenuation on all lines occurs inside a 20 dB directional coupler, which allows us to minimize the microwave power 
dissipated on the cold stage of the cryostat. The JPA is a reflection amplifier; a signal incident on the strongly coupled port of the JPA 
is reflected and amplified. A cryogenic circulator is used to separate the incident and reflected waves, defining the input and output 
ports of the JPA. The other circulators are used to isolate the cavity from the noise emitted from the amplifier's input. 



where n' add is the total added noise of the measurement in 
units of equivalent number of photons. For an ideal mea- 
surement (i.e. for a quantum-limited measurement of both 
quadratures of the light field), n' add = 1/2. 



Before the onset of normal-mode splitting, one can di- 
rectly relate the measured microwave power spectrum S to 
the displacement spectral density S x . Assuming A = 0, 



n c <C n m and g,S -C k, 



<? 1 k r n T 

— = - + n ' + 4/3^r (S6) 

hw 2 add k (r m + rf + w 



1 2(3G 2 n d k cx 

2 +n add + K K ' 



(S7) 



where V = 4g 2 /K is the optomechanical damping rate. 
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MICROWAVE MEASUREMENT AND CALIBRATION 

The detailed circuit diagram for our measurements is 
shown in Fig. S2. In order to calibrate the value of go = 
Gx zp for this device, we applied a microwave drive op- 
timally red-detuned (A = 0) and measured the thermal 
noise spectrum of the mechanical oscillator as a function 
of cryostat temperature. Here we restricted rid ~ 3 in or- 
der to ensure that radiation pressure effects are negligible. 
With the value of g now determined, we increase the drive 
amplitude and measure the thermal noise spectrum at each 
drive power. The noise spectra are recorded and averaged 
with commercial FFT spectrum analyser. Each spectrum 
is typically an average of 500 traces with a measurement 
time of 0.5 s per trace. The cavity response is then mea- 
sured with a weak probe tone with a vector network anal- 
yser to determine precise cavity parameters at each mi- 
crowave drive power, including the precise detuning and 
K. For larger microwave drive powers where the cavity 
spectrum exhibits optomechanically induced transparency 
effects [7? ? ], this spectrum also serves as a direct mea- 
sure of g. Finally, using additional calibration tones, each 
noise spectrum is calibrated in units of absolute microwave 
noise quanta and fit with Eq. 5 to determine the occupancy 
of both the cavity and mechanical modes. 

For our measurements, we infer that our entire measure- 
ment chain has an effective added noise of n' add = 2.1. 
This value is consistent with the independently measured 
value for the added noise of the JPA (n a( jd = 0.8) and the 
2.5 dB of loss between the output of the cavity and the JPA 
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FIG. S3. Measured transmission of filter cavity. A tunable 
resonant cavity was implemented at room temperature in order to 
suppress noise ~ 10 MHz above the drive frequency. As shown 
here, this cavity reduces the noise at the cavity frequency by more 
than 40 dB, ensuring that the phase or amplitude noise of the 
generator is not responsible for the finite occupancy of the cavity 
at large drive power. 

In order to ensure that the amplitude or phase noise of 
the signal generator was not responsible for the finite oc- 
cupancy of the cavity at high drive power, we designed and 



built a custom filter cavity ]4). As shown in Fig. S3, when 
the filter cavity is tuned to precisely the frequencies of our 
circuit, it provides an addition 40 dB of noise suppression 
at the cavity resonance frequency. The phase and ampli- 
tude noise of our signal generator alone are specified by the 
manufacturer to be less than -150 dBc at Fourier frequen- 
cies 10 MHz away from the drive. With the addition of fil- 
ter cavity, we lower this noise to well below the shot-noise 
level of our microwave drive. Furthermore, even without 
the filter cavity, we could not resolve an appreciable differ- 
ence in the cavity occupation. Thus, while we do not know 
the precise mechanism for this occupancy, we conclusively 
determine the generator noise is not the cause. 



INFERRING CAVITY PARAMETER AND NUMBER OF 
DRIVE PHOTONS 

The measured microwave cavity parameters may be in- 
ferred from the transmitted power spectrum. The power at 
the output of the cavity P is related to the input power P t 
byQ 



P = P 

1 O 1 1 



n 2 + 4A ; 



(S8) 



where A = u d — uj c is the difference between the fre- 
quency of the drive uj^ and the cavity resonance frequency 
uj c . k is the total intensity decay rate of the cavity (full 
width at half maximum) with k = k + k cx . k is the 
coupling rate to the dissipative environment, and k cx is the 
coupling rate to the transmission line used to excite and 
monitor the cavity. 

The number of photons in the cavity due to a coherent 
input drive at detuning A may be calculated from the stored 
energy E in the cavity. 



n d 



E 2P[ k cx 
fiWd huJd k 2 + 4A 2 



(S9) 



For our circuit, k cx = 2tt x 133 kHz. Thus, when the 
drive is optimally detuned such that A = — O m , the input 
power required to excite the cavity with one photon is Pj ~ 

2hu: d n 2 jK cx « 50 fw. 



FUNDAMENTAL LIMITS OF SIDEBAND COOLING 

Equation 2 in the main text gives an expression for the 
final occupancy of a mechanical mode, assuming that the 
microwave drive is optimally detuned (A = — f2 m ). This 
expression is only the lowest order approximation in the 
small quantities g/£l m and K/T2 m . Up to second order, the 
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final occupancy is [[T(| 



n n 



+ n c 



T m Ag 2 + n 2 
K Ag 2 + KT ] 

Ag 2 + «r E 



1 + 



g 2 Ag 2 + «r m 
Q m Ag 2 + k 2 
8g 2 + k 2 Ag 2 + kT z 



sn 2 



Ag 2 



+ 



g 2 + k 2 



The last term represents the fundamental limit for sideband 
cooling and demonstrates the importance of the resolved- 
sideband regime. For our system, Q m S> k, g; and hence 
this last term only contributes negligibly to the final occu- 
pancy of the mechanical mode (< 10~ 4 quanta). 



Measurement imprecision and backaction 

Throughout the main text and this supplementary infor- 
mation, we use the "single- sided" convention for all spec- 
tral densities in which for any quantity A, the mean-square 
fluctuations are {A 2 } = / °° S a {uj)^. This yields the fa- 
miliar classical result that an oscillator coupled to a thermal 
bath of temperature T will experience a random force char- 
acterized by the force spectral density Sp = Ak B TmT m . 
More generally, 



S F = AhQ n 



T I 



(S10) 



where n m is the Bose-Einstein occupancy factor given by 



[exp(h£l m /k B T) - 1]~ 



Independent of any convention for defining the spec- 
tral density, the visibility of a thermal mechanical peak of 
given mechanical occupancy above the noise floor of the 
measurement represents a direct measure of the overall ef- 
ficiency of the detection. As shown in Fig. S4, ratio of 
the peak height to the white-noise background allows us 
to quantify the imprecision of the measurement in units of 



mechanical quanta ]9], n 
spection of Eq. S6 implies 



imp 



S™ p mn m T'J{Ah). In- 



n-y 



1 k Ag 2 + KT m ( 1 



4/3 k c 



Ag 2 



+ < 



add 



(Sii) 



Once the drive is strong enough that g 3> y / Kl\^), n imp 
no longer decreases with increasing drive. It is precisely 
because we are measuring with a detuned drive that also 
damps the mechanical motion, that nj mp asymptotically ap- 
proaches a constant value l3l [TT1 . For an ideal measure- 
ment OS = 1,k = k cx , and < dd = 1/2), n imp -> 1/4. 
Implicit in obtaining this optimal value for n' add and hence 
nimp is that all the photons exiting the cavity are measured. 
Any losses between the cavity and the detector can be mod- 
eled as a beam-splitter that only transmits a fraction ij of the 



photons to the detector and adds a fraction (1 — rj) of vac- 
uum noise. So the effective added noise n' add accounting 
for these losses becomes 



n 



add 



n add 
7] 



V 



(S12) 



Thus, shot-noise limited detection of the photons (n add = 
1/2) is a necessary, but not sufficient, condition for reach- 
ing the best possible level of precision. 

Quantum mechanics also requires that a continuous dis- 
placement measurement must necessarily impart a force 
back on the measured object. For an optimally detuned 
drive (A = 0) in the resolved-sideband regime, this back- 
action force spectral density S F a approaches a constant 
value as a function of increasing drive strength and asymp- 
totically approaches S F a = 2hD, ai mT' m . Again, express- 
ing the spectral density in units of mechanical quanta gives 
n ba = S£ a /(4/^ m mr m ) -> 1/2. 

Fundamentally, the Heisenberg limit does not restrict the 
imprecision S" up or the backaction alone, but rather it 
requires their product has a minimum value ll3l lT2l 



Qimp Qba 



imp 



> h. 



(S13) 



An ideal cavity optomechanical system can achieve this 
lower limit for a continuous measurement with a drive 
applied at the cavity resonance frequency. When con- 
sidering the case where the drive is instead applied de- 
tuned below the cavity resonance (A = 0), this product 
never reaches this lower limit J3j ITTj and is at minimum 

To estimate these quantities for our measurements, we 
can infer the total force spectral density experienced by 
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FIG. S4. Measurement imprecision in units of mechanical 
quanta. 
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our oscillator as Sf tal = 4Ml m mr^(n m + 1/2). As 
this total necessarily includes the backaction, we may make 
the most conservative assumption that it was solely due 
to backaction that our oscillator remained at finite occu- 
pancy. Hence, n ba < n m + 1/2. The low thermal oc- 
cupancies attained in this work allow us to place an upper 
bound on how large the backaction could possibly be, and 
hence quantify our measurement in terms of approach to 

the Heisenberg limit. Thus, \J S™ v S b p = 4fi^/n imp n ba < 

4/jy / n imp (n m + 1/2). At n d = 3 x 10 4 , we simulta- 
neously achieve n m = 0.36 and n imp = 1.9 (S^ otal = 
1.6 x irr 34 N 2 /Hz and S , ™ p = 1.7 x 10" 33 m 2 /Hz) 
yielding an upper limit on the measured product of back- 
action and imprecision of 5.1 h. As stated above, the best 
possible backaction-imprecision product is hs/2 when us- 
ing red-detuned excitation; thus our measurement is only a 
factor of 3.6 above this limit. It may also be noted that this 
factor would have been 1.8 except that our chosen geome- 
try losses half of the signal back to the input ([3 = 1/2). In 
future experiments, using a single-port geometry = 1) 
will improve this inefficiency. 
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